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Abstract. In this paper we derive an easily computed approximation of Rogers and Shi's 
lower bound for a local volatility jump-diffusion model and then use it to approximate European 
basket option values. If the local volatility function is time independent then there is a closed-form 
expression for the approximation. Numerical tests show that the lower bound approximation is 

' fast and accurate in comparison with the Monte Carlo method, the partial-exact approximation 

. method and the asymptotic expansion method. 
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1 Introduction 



It is in general difficult to value basket options due to lack of analytic characterization of the 
distribution of the underlying basket asset process. Research is mainly focused on developing 
fast and accurate approximation methods and finding tight lower and upper bounds for basket 
option values. In the Black-Scholes setting Curran (1994) and Rogers and Shi (1995) derive a 
lower bound for Asian options by conditioning random variables and Jensen's inequality. Deelstra 
et al. (2004) obtain the bounds for basket options with the comonotonicity approach. In affine 
^ Levy models one can derive lower bounds numerically for arithmetic Asian options based on 

, the characteristic function and the method developed in Duffie et al. (2000), see Albrecher et 
Tlj" I al. (2008). Deelstra et al. (2010) provide a good overview of the recent development in finding 
and computing the bounds. 

It is known that Rogers and Shi's lower bound is generally tight and is one of the most 
accurate approximations to basket option values. The lower bound can be calculated exactly 
CS| ' in the Black-Scholes framework. Xu and Zheng (2009) show that the lower bound can also be 
calculated exactly in a special jump-diffusion model with constant volatility and two types of 
Poisson jumps (systematic and idiosyncratic jumps). The usefulness of Rogers and Shi's lower 
^ ' bound depends crucially on one's ability of finding some highly correlated random variables to 
^ ■ the basket value and computing the conditional expectation exactly. It is difficult to extend 
Rogers and Shi's lower bound to more general models such as local volatility models due to 
lack of explicitly known distributions for models with non-constant volatilities, see Albrecher et 
al. (2008). To the best of our knowledge, Rogers and Shi's lower bound for models with local 
volatilities has not been discussed in the literature. 

In this paper we aim to find a good approximation to Rogers and Shi's lower bound for a 
local volatility jump-diffusion model and then use it to approximate the European basket option 
values. We first apply the second order asymptotic expansion (see Benhamou et al. (2009)) to 
approximate the basket asset value, then choose a normal variable and a Poisson variable as 
conditioning variables which are highly correlated to the basket asset value, and finally apply the 
conditional expectation results of multiple Wiener-Ito integrals (see Kunitomo and Takahashi 
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(2001)) to approximate Rogers and Shi's lower bound. The main contribution of the paper is the 
derivation of an approximation to Rogers and Shi's lower bound for local volatihty jump-diffusion 
models. We suggest an easily implemented algorithm to compute the lower bound approximation. 
If the local volatility function is time independent then there is a closed-form expression for the 
approximation. Numerical tests show that the lower bound approximation is fast and accurate in 
most cases in comparison with the Monte Carlo method, the partial-exact approximation method 
(Xu and Zheng (2009)) and the asymptotic expansion method (Xu and Zheng (2010)). 

The paper is organized as follows. Section 2 formulates the basket local volatility jump- 
diffusion model and explains some known methods in pricing European basket options. Section 
3 applies the second order asymptotic expansion to derive an approximation to Roger and Shi's 
lower bound which can be easily computed. Section 4 compares the numerical performance of the 
lower bound approximation with other methods in pricing basket options. Section 5 concludes. 

2 Model 

Assume {il.,P,J-,J-t) is a filtered risk- neutral probability space and J-j is the augmented natural 
filtration generated by correlated Brownian motions Wi,...,Wn with correlation matrix (pij) 
and a Poisson process N with intensity A. Assume Brownian motions Wi, . . . , Wn and Poisson 
process are independent to each other. Assume the portfolio is composed of n assets and the 
discounted asset prices Si, . . . ,Sn satisfy the stochastic differential equations 

dSiit) = a^{t, S^{t))dW^{t) + hA{t-)dM{t), (1) 

where /ij > — 1 are constant jump sizes and M{t) = N{t) — Xt is the compensated Poisson process. 
The discounted basket value at time T is given by 

n 

S{T) = ^wA{T), 
1=1 

where Wi are positive constant weights. The European basket call option price at time is given 
by 

C{T,K) = B[{S{T) - K)+] 

Almost all work in the literature on Asian or basket options valuation assume the underlying 
asset prices follow lognormal processes, which corresponds to ai{t,Si) = aiSi with ai being 
positive constant and hi = for all i in ([T]). Since Asian options are similar to basket options, 
we do not differentiate these two types of options, even though some techniques are originally 
developed for Asian options. We now review some well-known approaches in approximation and 
error bound estimation for the pure diffusion case. 

Levy (1992) approximates the basket value S{T) with a lognormal variable which has the 
same first two moments as those of S(T) and derives an approximate closed-form pricing formula 
for C{T, K). The result is very good when maturity T and volatilities cjj are relatively small. The 
performance deteriorates as T or cxj increases. Curran (1994) introduces the idea of conditioning 
variables. Assume A is a random variable which has strong correlation with ^(r) and satisfies that 
S{T) > K whenever h> d\ for some constant d\. The basket option price can be decomposed 
as 

E[(5(r) - K)+] = E[iS{T) - if)l[A>,^]] + E[{S{T) - if)+l[A<d^]]. 

Curran (1994) chooses A a normal variable (geometric average) and finds the closed- form expres- 
sion for the the first part and uses the lognormal variable and the conditional moment matching 



2 



to find the approximate value of the second part. Rogers and Shi (1995) use the conditioning 
variable A to derive the lower and upper bounds to the basket option price C{T,K) as 

LB = 'E[{'E[S{T)\A]- K)+] 

and 

UB = LB+^B [var(5(T)|A)l[A<,^]]5 E[l[A<rf^]]i 

These bounds may be computed analytically for specific conditioning variables A. Numerical 
tests show that the lower bound is normally very tight to the true (Monte Carlo) value. 

Xu and Zheng (2009) extend the partial exact approximation method of Curran (1994) to 
jump-diffusion processes ([T]) with local volatility functions ai{t,Si) = aiSi for all i. (The jump 
part of the model in Xu and Zheng (2009) is more general than ([1]) with both common and 
individual jumps.) Asset prices Si(T) at time T have closed form expressions 

Si{T) = Si{0) exp (-^^rfr - h,XT + aiWi{T) + N{T) ln(l + hi)^ , (2) 

where Si{0) are asset prices at time 0, for i = 1, . . . , n. The conditioning variable A is chosen to 
be 

A = mN(T) + aW{T), 

where N{T) is a Poisson variable with parameter AT, W{T) a standard normal variable, N{T) 
and W{T) are independent to each other, m = Yli=i ^^(l + ^i)> = Z^fj=i o.iO^jPij'^i'^jT, and 
Qi = WiSi{0) exp((— ^(jj? — hi\)T) for z = 1, . . . , n. We then have S{T) > K whenever K> K — c, 
where c = ^11=1 ^i- N(T) = k and W{T) = y it is easy to find 

n ,^ 

E\S{T)\k = mk^ay\=Y^ ai exp ( -{a^T - Rf) + Riy + k ln(l + h^ 

i=i ^ 

where Ri = ^ Sj=i O'jPij^i^j'^- The lower bound can be computed from the expression 

LB = Y,Pk {E[S{T)\A = mk + ay] - K)+d^y), 
k=o ■'-'^ 

where pk = P{N(T) = k) and <^ is the cumulative distribution function of a standard normal 
variable. Xu and Zheng (2009) suggest a partial exact approximation (PEA) to the basket option 
value C (T, K) by the following value 

3 

C\T, K) := B[{S{T) - K)+1[A>^_,]] + g.E[(E[5(T)| A] + - K)+1[^^k-c]1 

i=l 

where qi = 1/6, q2 = 2/3, = 1/6, and ai = — \/3eo, a2 = 0, as = \/3eo, eo is some constant 
depending on the conditioning variable A and the conditional second moment of ^(r) given 
A. The lower bound plays a dominant role in the approximation with a weight 2/3, the other 
two parts with a weight 1/6 each are the second moment adjustment to the lower bound. Xu 
and Zheng (2009) show that LB < C^{T,K) < UB, which provides the error bounds for the 
approximate basket option value. The PEA method is fast and accurate in comparison with some 
other well known numerical schemes for basket options, see Xu and Zheng (2009) for details. A 
key condition for the PEA method is that one must know the closed form expressions of asset 
prices Si{T) as in ([2]) for i = 1, . . . ,n. This is impossible for general local volatility functions. 

Xu and Zheng (2010) discuss the basket options pricing for asset model ([T]) with the same jump 
sizes hi = h for i = 1, . . . , n. Since the PEA method cannot be applied for general local volatility 
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functions o"j(t, Si), Xu and Zheng (2010) first reduce the dimensionahty of the problem by directly 
working on the portfolio asset value S{t) which satisfies the following stochastic volatility jump 
diffusion model 

dS{t) = V{t)dW{t) + hS{t-)dM{t) 

with the initial price S'(O) = Yll=i ^iSiiO), where is a standard Brownian motion, independent 
of M, and 

n 

^iif ■= X] WiWjai{t,Si{t))aj{t,Sj{t))pij. 

The European basket call option price C{T,K) at time satisfies a partial integral differential 
equation (PIDE) 

Ct{T, K) = XhKCKiT, K) + ia(r, KfCKK{T, K) + \{h + 1) {c{T, - C{T, K) 

with the initial condition C(0, i^) = (5'(0) — K)'^ and the local variance function 

a{T,Kf = B[V{Tf\S{T) = K]. 

The main difficulty is to compute the conditional expectation as asset prices Si{T) have no 
closed-form expressions. Following the first order asymptotic expansion methods of Kunitomo 
and Takahashi (2001) and Benhamou et al. (2009), Xu and Zheng (2010) show that the unknown 
local variance function a(T, K)"^ can be approximated by 

a{T,Kf^a{T) + h{T){K-Sm, (3) 

where a{T) = E"i=i WiWjPijPiPj, h{T) = E"j=i j^WiWjPijPiPj i^fCi + |^Cj) , pi = ai{T, S'i(O)), 

* = as|oy'^»(^''^»(0))' Ci = YJj=iWj Pijijg ai{t,Si{0))aj{t,Sj{0))dt) , and = EiLi^iC'i- 
One can find the basket option price by solving the PIDE with the implicit-explicit finite difference 
method. The method is less accurate than the PEA method when the local volatility function is 
the Black-Scholes type, but is more general and can handle general local volatility functions. The 
main limitations of the method are that the approximate variance function in ^ is not always 
positive, which can cause numerical errors in solving the PIDE, and that it requires the same 
jump size for all assets, which is unrealistic for a general asset portfolio. 



3 Lower Bound Approximation 

The asymptotic expansion for basket options pricing in a general diffusion model is discussed in 
Kunitomo and Takahashi (2001) in which the valuation of conditional expectations is a necessary 
step to obtain the characteristic function of the basket value. In this paper we use the same 
method to expand the parameterized asset price processes to the second order, see Benhamou et 
al. (2009), and apply the conditional expectation results of multiple Wiener-Ito integrals directly 
to approximate the Rogers and Shi's lower bound for the jump-diffusion model ([1]). For e G [0, 1] 
define 

dSfit) = eai{t, S'i{t))dWi{t) + eh,Sl{t-)dM{t) (4) 
with initial condition 5|(0) = S'i(O). Note that SjiT) = Si{T). Define 

(fc) W ._ aV.(t,5|(t)) 

' ^ ' agfc 1^=0 a'^^ — — Q^sff — 
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for A; = 0, 1, .... It is obvious that ^ (t) = Sf{t) = Si{0) and cif ^ (t) = ai{t, Si{0)) for for alU > 
and i = 1, . . . , n. In particular, if local volatility functions ai are time independent then functions 

(k) 

o") are all constants for i = 1, . . . ,n, which may simplify considerably some computations. An 
important case is the constant elasticity of variance (CEV) model, where ai{t,Si) = UiSf' for 
i = 1, . . . , n. To simplify the notations in subsequent discussions we denote by 

n 

af\t) = Y,w,af\t)<jf\t)p,,, k = 0,l. 

The second order asymptotic expansion around e = for S^{t) is 

St{T) « sf\T) + sl'\T)e + ^^f )(r)6l 
Expand (H]) to the second order, we have 

dSl^\t) = af\t)dWiit) + hiSi{0)dM{t), 

dSf\t) = 2af\t)Sf\t)dW^{t) + 2hiS^p{t-)dM{t), 

with the initial conditions 5j^^^(0) = sf'\^) = 0. Therefore, 

S^^\t) = -\hiSmt+ f af\s)dWi{s) + hiSmN{t) (5) 







Sf\t) = -2\h^ f sl'\s)ds + 2 f al'\s)si'\s)dWi{s) + 2h^ f sl'\s-)dN{s) (6) 
Jo Jo Jo 

for < t < T. Letting e = 1 we may approximate the basket value S{T) by 

S{T) ^ S^{T) := S{0) + S^^\t) + ^^-^, (7) 

where M^HT) := Y.7=i WiSi^\T) for j = 1, 2. 

Since there are no closed- form expressions for Si{T) and S{T), it is difficult to compute the 
lower bound. If we approximate S{T) by S^{T), defined in ([7]), we may be able to compute the 
conditional expectation E[5'^(T)|A] for some conditioning variable A and then to approximate 
the lower bound. We therefore propose that 

LB ^ LBA:=E[(E[S^{T)\A]- K)+]. (8) 

The next step is to choose the conditioning variable for the approximation. In the Black-Scholes 
framework, it is without exception to choose normal random variable as the conditioning vari- 
able, as the asset price is log-normally distributed. It is not clear what one should choose 
for general jump-diffusion models as there are no closed form solutions for underlying asset 
prices. With the insight from Xu and Zheng (2009) we choose A = (A^(r),A(r)) for the 
model, where A(r) = E]=iWj Jq af {t)dWj{t) is a normal variable with mean and vari- 
ance = J2i=i'^i lo ^i^\'^)d't, and N(T) is a Poisson variable with parameter AT. From ([8]) 



and d?]) we have 

LBA=J2Pk 5(0) +E[5«(T)|A = {k,vx)] +E[^^\A = {k,vx)] - k] 

k=o ■'-^ LV ^ J 



d^{x). 
(9) 
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Since ^(^^(r) = Y^i=lW^s\^\T), we only need to find E[5p'^(r)|A = {k,vx)] for j = 1,2 and 
i = 1, . . . ,n. We first calculate E[S'f ^(r)|A = {k,vx)]. 



nsf\T)\k = {k,vx)] 

- XhiSi{0)T + B[hiSiiO)N{T)\N{T) = k] + B 



af'\t)dW,{t)\A{T) = vx 



XhiSi{0)T + hiSi{0)k + 



T 



df\t)dt\ X 



The valuation of E[;^^2^|A = {k , vx)] is more involved. Since ^(T) is the sum of three 

terms in ([6|) we may find the conditional expectation of each term by substituting S^^\t) in ([5]) 
into the integrands and then computing three conditional expectations. We now derive them one 
by one. The first term can be written as 



T 



5f^(t)dt|A = {k,vx) 



A1+A2 + A3 



where 



^1 



-Xhi)B 



-XhiSiiO)t)dt\A = ik,vx) 







-S,{0)X'h^t 



2 ^2r7-i2 

T ft 



A, 



-A/i,)E 
-A/ii)E 



JO 
T 



af\s)dWi{s))dt\A = {k,vx) 



{T-t)af\t)dWi{t)\A{T) = vx 



-Xhi)'^(^j\T-t)af\t)dt^x 

-A/i,)E f {hiSi{0)N{t))dt\A = {k,vx) 
Jo 

-XS^{0)hf) [ B[N{t)\N{T) = k]dt 
-XS,{0)hj 



Here we have used the fact that {N(t)\N{T) = A;) is a binomial variable with k independent 0-1 
trials and probability j 
term can be written as 



trials and probability ^ = ^ of taking 1, which implies E[A^(t)| A^(T) = k] = y- The second 



E 



al'\t)sl'\t)dWi{t)\A = {k,vx) 



B1 + B2 + B3, 
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where 



Bi 



Bo 



B^ 



E 



af\t) {-XhA{0)t) dWi{t)\A = {k,vx) 



T 



i-Xhi)S,{0){ I M^\t)dt]^x 



E 



( / al'"{s)dWi{s) ) dWi{t)\A = {k,vx) 



Ml 



T / rt 



\J0 



a^'\s)ds]al'\t)dt]{x'-l) 



E 



4'\t) {hA{0)N{t))dWi{t)\A = {k,vx) 



j\f\t)^dWm^{T) = vx 



k 



tal^\t)dt^^x. 



The computation of B2 is discussed in Kunitomo and Takahashi (2001, Lemma A.l). The third 
term can be written as 



where 



sl'\t-)dNit)\A = {k, 



vx] 



Ci + C2 + C3, 



Ci = hiE 



XhiSi{0)?jdN{t) 



A = (k, vx) 



--SiiO)XhfTk 



Co 



hiE 



T 



\J0 
T / ft 



(jf\s)dWi{s) ]dN{t) 



A = {k, vx) 



hiE[j^ (^j\f\s)dW,{s)^^dt\^{T) = vx] 



k 



t- rj-] 

C3 = hiE 



T 



{T-t)af'\t)dt]^x 



T 



h]SME 



hiSi{0)N{t-)jdN{t) 

Af(T)-l 

E mT) = k 

1=0 



A = {k, vx) 



hfSM 



k^ -k 



Substituting the first and second order conditional expectations into ([9]) we get the lower 
bound approximation as 



°° r-OQ / \ + 

LBA = '^Pk icx^ + ai{k)x + ao{k)] d^{x), 



(10) 
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where 



n 1 " 

ao{k) = S{0) + {K- AT) ^ WihiSi{0) + - ^ «;i5i(0)/i2 {{k - \Tf -k)-c-K 

i=l 1=1 

Note that in computing the lower bound approximation (|1U|) there is no need to do numerical 
integration. This is because for every fixed k the integrand is the positive part of a quadratic 
function which contains only three cases: no root, one root, and two roots. Therefore the compu- 
tation of the lower bound approximation LB A is easy and fast. This is one key advantage over 
the other methods such as simulation or finite difference for the PIDE. 

4 Numerical Results 

In this section we conduct some numerical tests for European basket calls with underlying asset 
price processes satisfying ([T]) to test the performance of the lower bound approximation. The 
Monte Carlo simulation provides the benchmark results. The control variate technique is adopted 
to reduce the standard deviations. The following data are used in all numerical tests: the number 
of assets in the basket n = 4, the portfolio weights of each asset wi = 0.25 for i = 1, . . . ,n, the 
correlation coefficients of Brownian motions pij = 0.3 for i, j = 1, . . . ,n, the initial asset prices 
Si{Q) = 100 for i = 1, . . . , n, the exercise price K = 100, and the jump intensity A = 0.3. 

Tables [1] displays the European basket option values with four different methods: the Monte 
Carlo (MC), the partial exact approximation (PEA), the asymptotic expansion (AE) and the 
lower bound approximation (LBA). The local volatility function is a{t,S) = 0.25. We perform 
numerical tests for three constant jump sizes hi = h = —0.2212, —0.1175, —0.0606. It is clear that 
all three approximation methods perform well with relative errors less than 1%. We have done 
other tests with the same data except the volatility function being changed to a{t, S) = 0.55. 
The overall relative error is 0.6% for the PEA, 4% for the AE and 1.7% for the LBA. Tables [1] 
shows that the performance of the PEA is the best while the LBA performs better than the AE. 
The results suggest that the PEA is the best approximation method for the European basket call 
pricing when local volatility functions are of the Black-Scholes type. On the other hand, the AE 
and LBA are much more fiexible and can handle general local volatility functions. The LBA has 
the additional advantage of having closed-form solutions when the local volatility functions are 
time independent and can deal with different jump sizes to common jumps. The AE is slow in 
solving the PIDE with the finite difference method and requires the same jump size of all assets 
to common jumps. 

Table [2] displays the numerical results with the MC, AE and LBA and with different maturities 
(r = 1,3) and local volatility functions (c(5) = aS^ with a = 0.2,0.5 and /3 = 1,0.8,0.5). The 
jump size is hi = h = —0.2212 for all assets. The last row displays the average errors of the AE 
and LBA. Table [2] shows that the performance of the LBA is excellent with the average relative 
error 0.4%. The overall performance of the LBA is better than that of the AE, especially when 
the local volatility function is cr{S) = 0.55. Matlab is used for computation. The LBA only takes 
a few seconds for each case, much faster than the AE and MC. We have done other tests with the 
same data as in Table [2] except different intensity rates (A = and 1). The overall performance 
is essentially the same as that of Table [2j We can say with reasonable confidence that the LBA 
suggested in this paper works well. 



8 



Tables [3] displays the numerical results with the MC and LBA and with different maturities T 
and local volatility functions a{t, S) = aS^ . The jump sizes for four assets are hi = 0, /12 = 0.3, 
/13 = —0.3, and /i4 = 0, which implies that the jump event has no impact to asset prices 1 and 
4, positive impact to asset price 2 and negative impact to asset price 3. Since The PEA requires 
the Black-Scholes setting while the AE requires the same jump size for all assets, neither method 
can be used in this test. It is clear that the LBA performs well in comparison with the MC with 
the overall relative error less than 1%. 

5 Conclusion 

In this paper we discuss Rogers and Shi's lower bound approximation (LBA) to basket option 
values for local volatility jump-diffusion models. We expand parameterized asset prices to the 
second order with the asymptotic expansion method and obtain an easily implemented LBA. It 
turns out that if the local volatility function is time independent, such as the CEV model, then 
there is a closed-form expression for the approximation. We compare the numerical performance 
of the LBA with other methods using different parameters and show that the LBA is fast and 
accurate in most cases. 
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A 


h 


T 


MC (stdcv) 


PEA (crr%) 


AE (err%) 


LBA (err%) 


0.3 


-0.2212 


1 


7.35 (0.01) 


7.35 


(0.0) 


7.35 


(0.0) 


7.37 


(0.3) 






:-! 


12.93 (0.01) 


12.92 


(0.1) 


12.8.") 


(().()) 


12.8() 


(0.5) 




-0.1175 


1 


6.08 (0.01) 


6.08 


(0.0) 


6.07 


(0.2) 


6.09 


(0.2) 






3 


10.57 (0.01) 


10.56 


(0.1) 


10.49 


(0.8) 


10.57 


(0.0) 




-0.0606 


1 


5.66 (0.01) 


5.66 


(0.0) 


5.65 


(0.2) 


5.67 


(0.2) 






3 


9.83 (0.01) 


9.82 


(0.1) 


9.74 


(0.9) 


9.86 


(0.3) 


1 


-0.2212 


1 


10.78 (0.01) 


10.77 


(0.1) 


10.78 


(0.0) 


10.82 


(0.4) 






3 


18.64 (0.01) 


18.63 


(0.1) 


18.57 


(0.4) 


18.91 


(1.4) 




-0.1175 


1 


7.28 (0.01) 


7.28 


(0.0) 


7.28 


(0.0) 


7.31 


(0.4) 






3 


12.65 (0.01) 


12.64 


(0.1) 


12.58 


(0.6) 


12.68 


(0.2) 




-0.0606 


1 


6.02 (0.01) 


6.02 


(0.0) 


6.01 


(0.2) 


6.03 


(0.2) 






3 


10.45 (0.01) 


10.43 


(0.2) 


10.37 


(0.8) 


10.47 


(0.2) 


Average 






(0.1) 




(0.4) 




(0.4) 



Table 1: The comparison of European basket call option prices with the Monte Carlo (MC), 
the partial exact approximation (PEA), the asymptotic expansion (AE) and the lower bound 
approximation (LBA). The table displays the results with different jump intensities A, jump sizes 
hi = h, maturities T, and local volatility function ai{t,S) = 0.25. The numbers inside brackets 
in the MC columns are the standard deviations and those in the PEA, AE and LBA columns are 
the relative percentage errors in comparison with the MC results. 



T 


a 




MC (stdcv) AE (err%) LBA (err%) 


1 


0.2 
0.5 


1 


7.35 (0.01) 7.35 (0.0) 7.37 (0.3) 
14.71 (0.01) 14.42 (2.0) 14.87 (1.1) 


0.2 
0.5 


0,8 


5.31 (0.01) 5.33 (0.4) 5.31 (0.0) 
7.33 (0.01) 7.33 (0.0) 7.34 (0.1) 


0.2 
0.5 


0.5 


5.09 (0.01) 5.09 (0.0) 5.08 (0.2) 
5.11 (0.01) 5.12 (0.2) 5.11 (0.0) 


3 


0.2 
0.5 


1 


12.93 (0.01) 12.85 (0.6) 12.86 (0.5) 
25.69 (0.04) 24.14 (6.0) 26.16 (1.8) 


0.2 
0.5 


0.8 


9.61 (0.01) 9.64 (0.3) 9.63 (0.2) 
12.86 (0.01) 12.86 (0.0) 12.81 (0.4) 


0.2 
0.5 


0.5 


8.96 (0.01) 8.98 (0.2) 8.91 (0.6) 
9.18 (0.01) 9.21 (0.3) 9.18 (0.0) 


Average 


(0.7) (0.4) 



Table 2: The comparison of European basket call option prices with the MC, AE and LBA. The 
table displays results with different maturities T, local volatility functions ai{t,S) = aS^ and 
jump sizes hi = h = —0.2212. 
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T 


a 


/3 


MC (stdcv) LBA (crr%) 


1 


0.2 
0.5 


1 


5.53 (0.01) 5.52 (0.2) 
13.87 (0.01) 13.95 (0.6) 


0.2 
0.5 


0.8 


2.22 (0.01) 2.22 (0.0) 

5.50 (0.01) 5.49 (0.2) 


0.2 
0.5 


0.5 


0.63 (0.01) 0.63 (0.0) 
1.42 (0.01) 1.42 (0.0) 


3 


0.2 
0.5 


1 


9.68 (0.02) 9.66 (0.2) 
24.42 (0.06) 24.84 (1.7) 


0.2 
0.5 


0.8 


3.95 (0.01) 3.94 (0.3) 
9.57 (0.02) 9.59 (0.2) 


0.2 

0.5 


0.5 


1.37 (0.01) 1.37 (0.0) 
2.59 (0.01) 2.59 (0.0) 




(0.3) 



Table 3: The comparison of European basket call option prices with the MC and LBA. The table 
displays results with different maturities T, local volatility functions ai{t,S) = aS^, and jump 
sizes hi = 0, /i2 = 0.3, hs = —0.3, and /14 = 0. 
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